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Abstract 

In the model of a first order quark- gluon/hadron phase transition in which 
the hadronic phase is considered as vacuum bubbles growing in the quark- 
gluon background with chiral symmetry broken inside the bubble, we find the 
estimate for the length scale associated with inhomogeneities originated during 
the transition, 10 m ~ ^ ~ 40 m, being sufficient to produce significant effects 
on cosmological nucleosynthesis. 
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Recent lattice gauge simulations indicate that the pure SU{?>) gauge theory has a first 
order phase transition between the low temperature color confining hadronic phase 
and the high temperature deconfined plasma phase Because of the first order 
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nature of the phase transition, these two phases can coexist at the critical temperature 
with a distinct phase interface in between; this interface carries a positive amount of 
extra free energy, i.e. a surface tension (a). 

The surface tension is a fundamental parameter in the description of the time 
evolution of hot hadronic matter as it proceeds through the phase transition. In the 
scenario of a homogeneous nucleation p[ the probability for a bubble to nucleate is 



p(T) = exp 



In this model, the bubble at nucleation time has a radius large enough to be hydro- 
statically stable in the supercooled plasma, and to sustain a difference in pressure 
{p — p) between two phases. 

The question has been raised in the literature 0, Q whether baryon inhomo- 
geneities could originate during a first order quark-gluon/hadron phase transition that 
would alter the successful conventional calculations of cosmological nucleosynthesis 
0. An essential ingredient in these calculations is the surface tension for bubbles of 
the new hadronic phase. The length scale associated with the inhomogeneities has 
been given in ref. [^] [| (we will arrive at this relation below): 

i ~ (3.4 X 10^ m)(a/MeV3)^/2(Tc/MeV)-^3/^ (2) 

where Tc is the critical temperature for a first order quark-gluon/hadron phase tran- 
sition. Campbell et al. |0] have obtained an estimate for this surface tension in the 
framework of the effective low energy Lagrangian for broken chiral and scale invari- 
ance p in the hadronic phase, 

a ^ (70 MeV)^ (3) 

which is in good agreement with the estimate made in ref. Recently the surface 
tension has been calculated in quenched QCD by lattice numerical methods [|, [11 



as well as analytically in the mean- field approximation [[12i p!3| . When performed on 



Nt = 2 lattices, the numerical studies have produced the results 

a/T^ = 0.24(6) (4) 



and |0| 



a/T^ = 0.12(2), (5) 



while for Nt = A lattice 



a/T-" = 0.027(4). (6) 

With a given by (3)- (6) and ~ 100 K, the formula (2) gives the estimate i ^ 
0.5 — 1 m, which is hardly of interest for nucleosynthesis since significant effects on 

^Eq. (2) takes into account the correction given in ref. H. 
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nucleosynthesis require £ > 10 m and probably ^ > 100 m.[| Because of a strong 
dependence of Tc, £ will be smaller for larger Tc, e.g., for Tc ~ 150 K, Eq. (2) with a 
given by (3) yields £ ~ 0.1 m. 
Bhattacharya et al. 



13] have found the expression for the interface tension be- 



tween distinct Z{N) vacua above the deconfinement transition for a pure SU{N) 
gauge theory (A^ = 2, 3), 



3V3N g 

where g is the gauge coupling constant. For a typical value of g 1 



0.1 



2 (corresponding 

0.3) at temperature of the order of a typical deconfinement 



to as = g'^/4:7T 

one, T ~ 150 MeV |T5[, the formula (2) yields the estimate 



£ ~ (15 - 30) m, (8) 

which is considerably larger than that obtained from Eqs. (3)-(6) above. 

In this letter we show that the estimate for i following from the values of the 
surface tension obtained in refs. [|], ^ ^ |Tl|, i ~ 0.5 — 1 m, should be raised 
by one-two order of magnitude, in agreement with Eq. (8), if one considers the 
new hadronic phase as vacuum bubbles growing in the quark-gluon environment, 
with chiral symmetry broken (i.e., nonzero value of the quark condensate) inside the 
bubble. 

At temperatures well below T^, hadronic matter consists of the lightest hadrons, 
the pions. At T near Tc, the vrvr interaction becomes strong, the vrvr amplitudes at 
the relevant energies become so large that they allow for some bound states, the 



resonances 



The method for taking into account such resonance interaction was 



suggested by Belenky and Landau |T3] as considering the unstable particles on an 
equal footing with the stable ones in the thermodynamic quantities, by means of a 



resonance spectrum. Such a spectrum in both the statistical bootstrap model ||T^, ^ 
and the dual resonance model 



takes on the form 



p(m) ~ m" 6"^/^°, 



(9) 



where a and To are constants. This treatment of a hadronic resonance gas leads to 
a singularity in the thermodynamic functions at T = To [|1^, |19[ and, in particular. 



to an infinite number of the effective degrees of freedom in the hadronic phase, thus 
making a transition to the quark-gluon phase impossible. To cope with this difficulty, 
it is normally suggested that the hadrons are strongly interacting particles with finite 
range interactions among them which increase with increasing density. These inter- 
actions cannot be ignored in discussing the thermodynamics of hadronic matter; in 



fact, it has been shown [21] that the neglect of such interactions leads to unphysical 



^As shown in ref. 
universe. 



m is needed for the most interesting models of a non-homogeneous 
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behavior for a deconfinement phase transition mentioned above. (The mass spectrum 
(9) originates from the clustering of hadrons such a clustering corresponds to an 



effective attractive interaction and would not allow for a transition to the quark-gluon 
phase since for this transition a repulsive interaction at small distances is necessary.) 
In ref. the hadronic resonance spectrum (9) was modified to include an effective 
repulsive interaction by means of i) a hard-core potential, ii) an excluded volume (Van 
der Waals approach). It has been shown that in both cases the number of the effective 
degrees of freedom in the hadronic phase is now finite (namely, 6.5-7, as read off from 
Figs. 6,7 of ref. [^), and a normal (first order) transition to the quark-gluon phase 
becomes possible. For our present purposes we shall restrict ourselves to the pions 
alone (i.e., to the 3 effective degrees of freedom), and think of the hadronic phase as 
that of the pions which are brought to the boiling point at Tc with the energy density 
p^(T,) = 37rV90 

The usual scheme of a first order phase transition is as follows. With the increase 
of energy deposition (e.g., in relativistic heavy ion collisions), bubbles of the quark- 
gluon plasma "vapor" are formed until at p ~ g/3 PniTc) the system is entirely 
vapor [g being the total number of degrees of freedom in the quark-gluon phase, 
which is 37 for two-flavor plasma and 47.5 for three-flavor one). The quark-gluon 
to hadron transition in the early universe proceeds in a slightly different way: when 
the universe is about 10 psec old, its temperature drops down to Tc ~ 150 MeV. 
However, the phase transition does not start immediately but the interface tension 
causes supercooling by the amount of AT/Tc ~ 0.02(ct/T^^)^/^ |2^. After the period 



of supercooling the hadron bubbles start nucleating; this process gives rise to shock 
waves which reheat the universe back to T^, thus hindering the further formation 
of the hadron bubbles p2| (in a different scenario, bubble formation and growth is 



sufficient to begin reheating the system due to the release of latent heat [^). The 
transition continues because of the growth of previously created bubbles, until the 
system gets entirely hadronized. The average distance between the bubbles turns out 
to be i 

i ~ 10"^ Rnia/T^f'^ m, (10) 
where Rh is the horizon size at the transition time, 

3 ■ 10^ m 

In view of (10), (11), one obtains the estimate (2) for the scale of inhomogeneities in 
the hadron distribution which may have affected cosmological nucleosynthesis. 

The transition temperature is calculated by equating pressures. If the quark-gluon 
plasma is treated within the bag model, then 

V. = 3^T,^ = PQCP = - B (12) 
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{B being the bag constant, B ~ (245 - 255 MeV)"^ [g4[), giving the value of ~ 170 
MeV for two-flavor plasma and ~ 160 MeV for three-flavor one. 

In conventional bag model the difference between the quark condensate values 
on different sides of a bag surface is normally disregarded. To estimate the effect of 
the quark condensate discontinuity on the properties of hadrons, we treat the letter 
within the bag model and use the energy-momentum tensor T^'^ of the quark fields 
ipjlx) with massless quarks inside the bag. In zero order approximation of a power 
expansion in the QCD coupling constant as the following expression holds: 

= ^ E (rfi'd'^rf - {d''rf)ri^'}) , (13) 

where a = 1,2,3 is SU(3)c color index, and f = u,d, s is SU(3)/ flavor index. The 
linear boundary condition]^ |^5| 

lYn^^rfix) = rf{x) (14) 

on the bag surface with external normal corresponds to the relation tp^ipj = 
on the interior side of the bag surface. On the exterior side i/jjip'j ^ because the 



corresponding vacuum condensates have nonzero values [26]. In this case the quark 



field contributions to the bag energy density pq = contain not only conventional 
part corresponding to the quark kinetic energy but also contain the usually ignored 
surface part. The surface part arises from the contribution of the discontinuity of the 
quark condensate values on both sides of a bag surface [ p7[| : 



Es = -lsY.{rfrf) = (^vacS. (15) 

/ 

Here S' is a surface area of the bag, and 

a^ac^-^Eirfrj). (16) 



Typical "empirical" values discussed in the literature are [^, ^ 

(CO ^ mr,) ^ -(240 ± 25 MeV)^ (17) 
(ipsips) being of similar magnitude. The commonly adopted value of the quark con- 



densate for calculations within the framework of QCD sum rules is p9 

{iIj^i/j") = -0.0241 GeV^ ^ -(289 MeV)l (18) 

■^In the chiral bag model, the boundary condition reads ijf^nf^iljj{x) — e^'^'''^''^^ 'i})'j{x), where 
h ~ 1,2,3 is the isospin index, and couples the quark fields on one side of the bag boundary with 
the pion field on the other side in a highly nonlinear manner. We do not discuss this model here. 
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Sometimes people consider even higher values of the quark condensate, e.g. 

(CO = ir.r,) = -(287 MeV)^ 

(V'.>:) = -(306 MeVf, (19) 

~ -0.032 GeV^ ^ -(315 MeVf. (20) 



or |3T1 



We shall take 



(CO ^ (CO ^ (CO ^ -(270 - 280 MeV)^ (21) 



which lies somehow between the values provided by (17)-(20). Then the value of ayac 
calculated from Eq. (16) is 

a^ac ^ (245 - 255 MeV)l (22) 

This surface tension is strong in comparison with the vacuum pressure (the bag 
constant) p4 



q (T 1 

B = -Pvac = -{-^G^^C^n - 4 E^/(CC) ^ (245 - 255 MeV)^ (23) 

Here is a gluon field stress tensor and is a typical bare (current) quark mass. 

Now consider the hadronic phase arising as a result of a first order quark-gluon/had- 
ron transition in the form of vacuum bubbles growing in the bulk quark-gluon plasma. 
Such a vacuum bubble is a sperical surface with quarks and gluons outside it obey- 
ing the bag boundary conditions^ (14). Inside the bubble we have nonzero quark 
{■ilj'^{x)ipf{x)) and gluon {G°^^G°'^'') condensates. Therefore, one can apply arguments 
similar to the previous ones and obtain the same expression for the interface surface 
tension as above (Eq. (16)). Thus, the value of a given by Eq. (22) should be used as 
the realistic value for the bubble surface tension in the expression for the length scale 
associated with the inhomogeneities, Eq. (2). The exect value for the critical tem- 
perature of a quark-gluon/hadron phase transition is not known. In the framework of 
chiral perturbation theory, Gerber and Leutwyler have calculated the tempera- 
ture of a chiral symmetry restoration transition (to be associated with the hadron to 
quark-gluon one) Tc ^ 170 MeV in the SU(2)j chiral limit. The value of the critical 
temperature provided by lattice gauge simulations is currently Tc ~ 140 MeV. 



We, therefore, consider the value of the critical temperature of a quark-gluon/hadron 
transition to be in the temperature interval 



140 MeV ~ Tc ~ 170 MeV, (24) 



''The estimate that we obtained, a'^/'^ ~ B^^'^, was also suggested in ref. 
^The bag boundary conditions for gluons are discussed in ref. [p3[. 



6 



as granted by both chiral perturbation theory and the most recent lattice simulations 
(note that both values of Tc estimated above within the bag model lie in this temper- 
ature range). For the temperature interval (24) and the value of the bubble tension 
given by (22), Eq. (2) gives 

10 m ~ £ ~ 40 m, (25) 

which is one-two order of magnitude as much as the estimate following from the values 
of a found in refs. ^ |10|, 0. In view of the aforementioned estimate i > 10 



[or 100) m, this value of the length scale should be considered as being sufficient to 



produce significant effects on nucleosynthesis. Note that the estimates of refs. [0, ^ 
lie in the interval (25). 

By further adjustment of the values of cr and Tc the estimate (25) may be raised 
or diminished, respectively. For example, for Tc ~ 100 MeV and cr given by (22) one 
obtains £ ~ 200 - 250 m. 



The question has been raised in ref. whether inhomogeneities originating dur- 



ing a first order quark-gluon/hadron transition can produce cosmo logical beryllium 
or boron. In contrast to simple inhomogeneities like those with i ^ 0.5 — 1 m, inho- 
mogeneities associated with the length scale (25) are likely to produce cosmological 
Be or B. If this is really the case, it will have great implications for Big Bang nu- 
cleosynthesis. Cosmological Be or B are a signature for significant density variations 
p5| . Such variations could lead to planetary mass black holes [^] or quark nuggets 



38| that could serve as the dark matter of the universe. 



As remarked in |l|], another interesting point which may lead to a non-standard 
scenario of nucleosynthesis is that there may be a baryon number contrast in the 
hadron and the quark-gluon plasma phases. The thermodynamic advantage to place 
most of the baryon number into the quark-gluon phase was analyzed in ref. fl^ 



A large penalty in the free energy is paid when a unit of baryon number is placed 
in the hadron phase because of the large mass of the baryons, as compared to the 
phase of massless quarks for which there is no mass penalty in the free energy. This 
may not be the case if an effective baryon mass is temperature dependent and drops 
significanfty as T gets closer to Tc, as suggested in a series of papers by Brown and 
Rho [^. It is most probably that the estimate of ref. |T^ for the baryon number 



contrast in two phases, R ~ O(IOO), should be diminished by one order of magnitude. 
Even in this case, the baryon number contrast in two phases will still remain large. 

Two phenomena, the large inhomogeneities scale in the hadron distribution and 
the large baryon number contrast in the hadron and the quark-gluon plasma phases 
should lead to sizable effects on the standard scenario of cosmological nucleosynthesis. 
With the observations made in this letter, conventional calculations of cosmological 
nucleosynthesis should be amended by nontrivial details. New work on this and 
related subjects is in progress. 
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